Preprint submitted to 

Acta Phys. Hung. A 22/1 (2005) 000-000 



HEAVY ION 
PHYSICS 



Nuclear kr in d+Au Collisions from Multiparticle Jet 
Reconstruction at STAR 

Thomas Henry^ for the STAR Collaboration" 

* Texas A&M University, 
College Station Texas, USA 



Received 31 October 2005 

Abstract. This paper presents the most recent nuclear kx measurements from 
STAR derived from multiparticle jet reconstruction of d+Au and p+p collisions 
at t/s = 200 GeV. Since jets reconstructed from multiple particles are relatively 
free of fragmentation biases, nuclear kx can be measured with greater certainty 
in this way than with traditional di-hadron correlations. Multi-particle jet 
reconstruction can also be used for a direct measurement of the fragmentation 
function. 
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1. Introduction 

Many analyses of Au+Au collisions at RHIC probe the dense medium using jets by 
studying medium modifications of the jet yields and jet shapes [^. The results from 
d+Au and p+p collisions provide the baseline for medium modification studies in 
Au+Au collisions. For these simpler systems, jets can be examined in great detail 
via full jet reconstruction. 

This study presents the most recent d+Au nuclear kx measurements from full 
jet reconstruction at RHIC based upon year 2003 STAR data. The transverse 
and longitudinal jet shapes characterized by the jx, and z distributions from p+p 
collisions are also presented, kx measures the transverse momentum of a parton 
within a nucleon or nucleus, and its growth from p+p to nuclear collisions is believed 
to contribute to the Cronin effect 
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2. Experiment 

The STAR detector is well suited for investigating jet production at RHIC, due to 
the complete <f> coverage and large r\ coverage of the TPC and EMC. The STAR 
Barrel EMC data used in this study provide neutral energy measurements including 
7T° decay photons. During the 2003 RHIC run, only half of the full barrel EMC 
was installed. This first half of the detector consists of 2400 towers of 0.05 <g> 0.05 
in t] (g> 4> each, leading to full azimuthal coverage in < r\ < 1. The barrel EMC 
was read out in minimum bias events. It was also used to trigger on "high tower" 
events, where one of the towers was above a nominal energy threshold of Et > 2.5 
GeV. This "high tower" triggered event sample contains a much larger fraction of 
jets than the minimum bias event sample. 
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Fig. 1. Average difference between the simulated parton Et and the total recon- 
structed jet Et- The gray region outlines the 1-standard deviation spread of this 
average difference. 

This analysis uses two different algorithms that reconstruct a jet by capturing 
the spray of fragmenting particles in a geometric cone. One centers the cone on the 
most energetic hadrons in the event, while the other optimizes the cone direction to 
maximize the included energy. Only the hadron centered cone algorithm is robust 
enough at high multiplicity to be used for reconstructing jets in d+Au collisions. 

Measurement of jet energy requires corrections for charged particle energy de- 
position in the EMC, for the finite efficiencies of the TPC and EMC, and for the 
unmeasured energy carried by long-lived neutral particles (n, Kl, •••)• It is also 
possible for the jet reconstruction algorithm to miss soft particles. 

A sample of pythia 1 events processed by the STAR detector simulator was 
analyzed in the same fashion as the data. Since the originating parameters for each 
parton were recorded in the monte-carlo pythia sample, the energy scale response 
of the STAR detector could be quantified. Figure 1 shows the total energy scale 
correction for jet radii of 0.7 and 0.5 as a function of E T '(Uncorrected)- From here 
on, Et will exclusively refer to reconstructed transverse jet energy, corrected for 
detector energy response using the curves in Figure 1. 

1 PYTHIA version 6.203 and 6.205. 
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3. Inclusive Jet Studies 



j —A . 



>- 
a. 



STAR Preliminary 

Cone Radius = .5 
\|s = 200 GeV 



A p T >0.6 GeV/c 
▼ p T >2.0 GeV/c 



n 




j T (GeV/c) 

Fig. 2. The left panel shows jet jt, the transverse component of the hadron 
momentum {j>h) relative to the jet axis, for charged hadrons. The upper (lower) 
data and grey PYTHIA monte-carlo histograms refer to hadron ph > 0.6GeV/c (ph > 
2.0GeV/c). For p h > 2.0 GeV/c, RMS j T = 486 ± 36 MeV. The right panel shows 
Ph/ Et for charged hadrons. The black data points and monte-carlo histogram are 
for 6.0 < E T < 7.0GeV. The grey data points are for 7.0 < E T < 9.0GeV. The 
upper range PYTHIA monte-carlo curve is statistically consistent with the lower, and 
is not shown. 



Jets in this study are required to pass the following cuts: number charged 
hadrons > 1, Cone Radius=0.5, g" h e " tra ^ < -65, and jet thrust axis > 40(70)cm 
from the inside(outside) of the BEMC (roughly .2 < r\ < .65 depending on z- Vertex). 
When analyzing jet shape, the perpendicular component of the hadron momentum 
(j't), and the fragmentation fraction z, are generally used. Figure 2 shows jet j+ 
and ph/Ex distributions using these cuts. Jets with high hadron count are biased 
toward higher detector energy response. The p^/Et distribution, since it is a per 
hadron distribution, is biased toward jets with high hadron count. This means 
Ph/ET is not directly comparable to the theoretical fragmentation function Dh(z), 
due to the bias introduced by the high hadron count jets. However, the substantial 
agreement between pythia and the data demonstrates that the phenomenological 
parameterizations embodied by pythia continue to be valid at these jet energies. 

The di-jet opening angle is another basic observable obtained from jet events. 
For di-jets, A</> = n in leading order QCD. Gluon radiation broadens A<fi, which is 
measured by the per-parton intrinsic kx = Et sin A<fi. When reconstructing di-jets, 
one "trigger" jet with observed Et > 7 GeV is reconstructed from both neutral and 
charged hadrons using the high tower which triggered the event. The other "away" 
jet direction is reconstructed using charged particles only, allowing the "away" jet 
to range —0.5 < r\ < 0.5. The "away" jet is not used to estimate jet energy scale. 

Figure 3 shows the di-jet fc-r distribution comparison for p+p and d+Au colli- 
sions. The p+p collisions have been supplemented with d+Au minbias background 
such that the multiplicity distributions of the supplemented p+p events matches 
that of the d+Au "high tower" sample. Taking cr^ T(ofcs) = a 2 kT{pp) + <J 2 kT{nucl) , the 
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\J a k T (nuci) ~ I'^l 0-35±o's7 GeV/c in d+Au collisions. The major systematic 
uncertainties are the jet energy scale in p+p and d+Au, the detector resolution, fit 
uncertainties, and d+Au background multiplicity. The fit uncertainties are largest. 
The systematic error on the d+Au cx^ T r n ucl) due to the detector kr resolution is 
small. In contrast, it is potentially large for the absolute measurement of <Jk T {pp)i 
and this effect is currently under study. 
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Fig. 3. kx distributions at y/s = 200 GeV. d+Au kr is shown in black, and p+p 
supplemented with d+Au background is shown in grey. Here kx has been corrected 
for energy but not thrust axis resolution. 

4. Conclusion 

Jets have been reconstructed in both p+p and d+Au collisions. The p+p jets have 
been used to measure jj>, z and kr- The inclusive jet fa and z results compare well 
with the PYTHIA monte-carlo simulation studies. The d+Au &T„„ ci value, though 
having large uncertainty, appears to be smaller than previous experiments [0] at 
much lower energy. The (j+) and y/ {kj) results from full jet reconstruction agree 
well with similar measurements using di-hadron correlations 
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